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£j , 

C/j , We report on the solar diurnal variation of the galactic cosmic-ray intensity observed by the 

Tibet III air shower array during the period from 1999 to 2003. In the higher-energy event samples 

k/ ' (12 TeV and 6.2 TeV), the variations are fairly consistent with the Compton-Getting anisotropy due 

\ to the terrestrial orbital motion around the sun, while the variation in the lower-energy event sample 

(4.0 TeV) is inconsistent with this anisotropy. This suggests an additional anisotropy superposed at 
the multi-TeV energies, e.g. the solar modulation effect. This is the highest-precision measurement 
of the Compton-Getting anisotropy ever made. 

PACS numbers: 96.40.Kk, 96.40.Pq, 96.50.Bh 



I. INTRODUCTION 

The galactic anisotropy of the cosmic-ray intensity is 
expected to carry information about the origin and the 
propagation mechanism of the galactic cosmic rays, as it 
reflects the magnetic field in space through which cos- 
mic rays have traveled. The anisotropy can be observed 
as the daily variation of cosmic-ray intensity recorded 
by the ground based detector in the sidereal time (side- 
real daily variation). The amplitude of the sidereal daily 
variation so far reported is as small as 0.1% or less, while 
the amplitude of the temporal variation of the air-shower 
event trigger rate amounts to a few% at the experimental 



site. Therefore, we have to eliminate the temporal varia- 
tion which is due mostly to the atmospheric pressure and 
temperature effects. 

We still have no reliable theoretical constraints for the 
sidereal daily variation due to the galactic anisotropy. 
This also makes it difficult to evaluate the systematic er- 
ror contained in the observed variation. Accordingly, we 
first analyze the daily variation in the solar time, which is 
expected from the so-called the Compton-Getting (C-G) 
anisotropy p| due to the earth's orbital motion around 
the sun. As this effect is predicted based on a reliable the- 
ory, the positive observation of this effect will assure the 
reliabilities of both the measurement and analysis. When 
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a cosmic-ray detector on the earth moves with respect to 
the rest frame of the cosmic-ray plasma, the fractional 
intensity enhancement due to the C-G anisotropy is ex- 
pressed, as Eq. 

^ = (7 + 2)^ cos8, (1) 

with I denoting the cosmic-ray intensity, 7 the power- 
law index of the cosmic-ray energy spectrum, v/c the 
ratio of the detector's velocity to the speed of light, and 
9 the angle between the arrival direction of cosmic rays 
and the direction of motion of detector Q • The vertical 
viewing detector on the earth moving along the circular 
orbit around the sun scans various directions in space 
as the earth spins and records the C-G anisotropy as 
the solar diurnal variation in its counting rate with a 
maximum at 6:00 hours in the local solar time. The 
amplitude of this variation is calculated to be as small 
as 0.05% or less, depending on the geographic latitude of 
the experimental site. 

The solar diurnal variation also can be caused by the 
anisotropy due to the solar modulation of galactic cos- 
mic rays in the heliosphere ||. Yasue et al. £j re- 
ported from their underground muon experiments that 
the anisotropy due to the solar modulation superposed 
on the C-G anisotropy extended up to several 100 GeV 
(see also @). 

Above 10 TeV energies, there have been only a few 
positive observational results reported on the solar di- 
urnal variation due to the C-G anisotropy @ @) be- 
cause of very small amplitude of the variation. Cutler and 
Groom reported in 1986 the first clear signature of the 
C-G anisotropy for multi-TeV cosmic rays (~ 1.5 TeV) 
in the solar diurnal variation . The variation reported 
by them was in reasonable agreement with the sinusoidal 
curve expected from the C-G anisotropy, while the max- 
imum phase of the curve deviated from 6:00 hours by +2 
hours at 2 a significance. They attributed the deviation 
to the meteorological effect on the underground muon 
intensity. 

Therefore, the solar diurnal variation is so far pre- 
sumed to be free from the solar modulation effect at TeV 
energies. We will report on the solar diurnal variation 
due to the C-G anisotropy observed by the Tibet III air 
shower array. The analysis of the sidereal anisotropy will 
be published elsewhere. 

II. EXPERIMENT 

The Tibet air shower experiment has been successfully 
operated at Yangbajing (90.522° E, 30.102° N, 4300 m 
above sea level) in Tibet, China since 1990. The array 
constructed first in 1990 was graduall y u pgraded by in- 
creasing the number of counters 

M U3 111 0, 

and 

then the Tibet III array, used in the present analysis, 
was completed in the late fall of 1999. This array con- 
sists of 533 scintillation counters of 0.5 m 2 each placed on 



a 7.5 m square grid with an enclosed area of 22,050 m 2 
and each viewed by a fast-timing (FT) photo-multiplier 
tube. A 0.5 cm thick lead plate is put on the top of 
each counter in order to increase the array sensitivity by 
converting 7 rays into electron-positron pairs. 

An event trigger signal is issued when any fourfold co- 
incidence occurs in the FT counters recording more than 
0.6 particles, resulting in the trigger rate of about 680 Hz 
at a few- TeV threshold energy. We collected 5.4 x 10 10 
events by the Tibet III array during 918 live days from 
November, 1999 to November, 2003. After some simple 
data selections (software trigger condition of any four- 
fold coincidence in the FT counters recording more than 
0.8 particles in charge, zenith angle of arrival direction 
< 45° , air shower core position located in the array, etc), 
3.0 x 10 10 events remain for further analysis. 

The pointing accuracy (0.02°) and angular resolution 
(0.9°) of the Tibet III array can be directly checked by 
monitoring the Moon's shadow in the cosmic-ray flux at 
multi-TeV energies . 

The performance of the Tibet III array is also exam- 
ined by means of a full Monte Carlo (MC) simulation 
in the energy range from 0.3 to 1000 TeV. We used the 
CORSIKA version 6.004 code and QGSJET model 
|l6| for the generation of air shower events and the EPICS 
UV7.24 code 01 f° r the detector simulation of shower 
particles with scintillation counters, respectively. Pri- 
mary cosmic-ray particles are sampled from the energy 
spectrum made by a compilation of direct observational 
data. The primary cosmic-ray energy is estimated by 
EppT which is the sum of the number of particles/m 2 for 
each FT counter. According to the result of the simula- 
tion, S/9ft = 100 corresponds approximately to 10 TeV 
primary cosmic-ray energy [l8| . 



III. ANALYSIS 

The selected air shower events are subsequently his- 
togrammed in hourly bins in the solar local time (365 cy- 
cle/year), according to event time, incident direction and 
air shower size of each event. To check the seasonal 
change in the daily variation, we obtained the histogram 
for each month and corrected it for the observation live 
time varying month to month. 

The daily and yearly event rates vary by ±2% and 
±5% 14], respectively, clue mostly to the meteorological 
effect. To eliminate these temporal variations and dis- 
cuss the daily variation with very small amplitude (0.05% 
or less), we adopt the following East— West subtraction 
method. We first obtain the daily variation in the so- 
lar time for each of 'East' and 'West' (E- and W-) in- 
cident events referring to the geographical longitude of 
the incident direction of each event and then subtracts 
the variation in the W-incident events from that in the 
E- incident events. Dividing this difference by the hour 
angle separation between the mean E- and W-incident 
directions averaged over the E- and W-incident events, 
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we finally reach the "differential" variation at solar time 
frame. This method Q largely cancels out the meteoro- 
logical effect and the possible detector biases, which are 
expected to produce common variations for both the E- 
and W-incident events. 

A possible drawback of this method is that we obtain 
only the "differential" form of the physical variation and 
we have to reconstruct the physical variation by "inte- 
grating" the "differential" variation with respect to the 
local time. This also makes the direct error estimation 
difficult. Hereafter, we will make the statistical argu- 
ment on the basis of the "differential" variation u D(t)" 
and compare the physical variation ll R(t)" with model 
curves obtained by integrating the best-fitted curve to 
"£)(*)". 

The data are then divided into the 3 data samples ac- 
cording to the representative primary energy of 4.0, 6.2 
and 12 TeV. Each representative energy is calculated as 
the mode value of the logarithmic energy of each event 
by the MC simulation. 

The expected event rate at our experimental site is 
calculated by considering the effective area of Tibet III 
and the C-G anisotropy given in Eq.Q]with 7 of 2.7. The 
expected ".D(t)" is then calculated by applying the E— W 
method and compared with the data. 



IV. RESULTS AND DISCUSSIONS 

Figure ^ shows the average solar daily variations ob- 
served by the Tibet III together with the sinusoidal 
curves best-fitted to the data. Note that the variations in 
Fig- ^a) are the "differential" variations ("D(i)"s) and 
the maximum phase in each panel is shifted earlier by 
6 hours (1/4 cycles) from the corresponding actual daily 
variation. The amplitudes in this panel are also 7r/12 
times as small as those of the actual variations. 

The x 2- fitting results for the 3 data samples, assum- 
ing a sinusoidal curve, are summarized in Table [I] The 
variations of the higher-energy data samples (12 TeV and 
6.2 TeV) are consistent with those expected from the C-G 
anisotropy, while the variation of the lower-energy data 
sample (4.0 TeV) statistically deviates from the expected 
C-G curve at 5.3 a significance in phase and at 8.3 a sig- 
nificance in amplitude, respectively. Compared with the 
expected C-G anisotropy, the x 2 value for the 4.0 TeV 
data sample is calculated to be 90.5/24 degree of freedom 
(d.o.f.), suggesting that the variation of this data sample 
is statistically inconsistent with the C-G anisotropy. It 
is noted here that the variation in the 6.2 TeV data sam- 
ple is consistent with the smooth interpolation between 
those in the 4.0 TeV and 12 TeV data samples as well as 
with the C-G anisotropy within statistics. 

A significant spurious variation in the solar time can 
be produced from the seasonal change of the sidereal 
daily variation due to the galactic anisotropy, as the av- 
erage variation in sidereal time is expected to be a few 
times larger than the solar daily variation due to the 
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FIG. 1: Average solar daily variations recorded in the 3 data 
samples with representative energies of 4.0, 6.2 and 12 TeV. 
The expected variations due to the C-G anisotropy are shown 
by broken lines, while the sinusoidal curves best-fitted to the 
data are shown by solid lines. From the top shown are the 
"differential" variations "D(t)"s (a), the physical variations 
"R(t)"s (b) and the differences between u R(t)" and the ex- 
pected variation (c), respectively. The error bars are statisti- 
cal. 
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TABLE I: Amplitudes A and phases <f> of the diurnal vari- 
ations of the \ 2 -fitted results assuming a sinusoidal curve, 
where the error bars are statistical in the data. Errors in the 
expected C-G anisotropy are negligible. 



C-G anisotropy [TfJ. Figure |21 shows "D(t)"s distribu- 
tions in the anti-sidereal (364 c/y) and extended-sidereal 
(367 c/y) time frames, which are both statistically in- 
significant. The insignificant variation (0.0021±0.0008%) 
in anti-sidereal time ensures that the seasonal change of 
the solar daily variation (365 c/y) is negligible. The in- 
significant variation (0.0015±0.0008%) in the extended- 
sidereal time, on the other hand, supports that the sea- 
sonal change of the sidereal daily variation (366 c/y) is 
also negligible. These results indicate that the spurious 
variation contained in the average variation is small. We 
thus estimate that the contamination due to the spurious 
variation might be less than 20% of the C-G anisotropy. 

The 4.0 TeV data sample in Fig. ^ is further divided 
into the 3 sub-samples according to the representative 
energies to check possible trigger threshold biases. We 
confirmed that each of these 3 sub-samples is statistically 
consistent with one another and the trigger threshold bi- 
ases are unlikely. 

To compare with the result by Cutler and Groom [jj, 
which is integrated above the muon threshold energy at 
128 GeV and is now the only one available at multi-TeV 
energies, we integrate the 4.0 TeV, 6.2 TeV and 12 TeV 
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FIG. 2: The "differential" variations "D(i)"s in the local 
anti-sidereal time (left) and extended-sidereal time (right). 
From the top shown are the average variations recorded in 
4.0 TeV (a), 6.2 TeV (b) and 12 TeV (c) data samples, re- 
spectively. The amplitudes of sinusoidal curves best-fitted to 
all data (3.0 xlO 10 events) are 0.0021±0.0008% in the anti- 
sidereal time and 0.0015±0.0008% in the extended-sidereal 
time, respectively. The error bars are statistical. 



by Cutler and Groom within statistics. Thus, the de- 
viation from the C-G anisotropy measured by our obser- 
vation in the 4.0 TeV data sample may be diluted by 
the energy integration. With our high statistics, we can 
afford to demonstrate the energy dependence of the so- 
lar diurnal anisotropy at multi-TeV energies for the first 
time. 

In conclusion, we clearly observe the C-G anisotropy 
in the 12 TeV and 6.2 TeV data samples. This is the 
most precise measurement of the anisotropy ever made. 
The measurement also suggests an additional anisotropy 
superposed on the C-G anisotropy in the 4.0 TeV data 
sample, for instance, the anisotropy due to the solar 
modulation effect. If this extra anisotropy is actually 
due to the solar modulation, the observed deviation 
from the C-G anisotropy may vary according to the solar 
activity changing in every 11 years. As the data used 
in the present analysis cover only the solar maximum 
period, we should examine this hypothesis by continuing 
the observation for a full solar activity cycle. The data 
during the next solar minimum around 2006 might be 
especially interesting. Furthermore, it may be also very 
useful to lower the energy threshold of the array down 
to sub- TeV region for a better understanding of the 
phenomenon. 



data samples and give a fit to it by assuming a sinu- 
soidal curve again. The amplitude normalized by the 
C-G anisotropy in our work becomes 1.28±0.08 (stat.) 
which is not inconsistent with 0.73±0.21 (stat.) obtained 
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